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We study the m-component vector spin glass in the limit m — > oo on a Bethe lattice. The 
cavity method allows for a solution of the model in a self-consistent field approximation and for 
a pcrturbative solution of the full problem near the phase transition. The low temperature phase 
of the model is analyzed numerically and a generalized Bose-Einstein condensation is found, as in 
the fully connected model. Scaling relations between four distinct zero-temperature exponents are 
found. 



I. INTRODUCTION 

The spin glass with m-component spins in the limit 
of large m is an interesting member of the spin glass 
family for various reasons. First, in the limit N — ► oo 
it is replica symmetric 1 . This sets it apart from high- 
dimensional and mean-field spin glasses with Ising, XY 
or Heisenberg spins but makes it more related to low- 
dimensional spin glasses, provided the replica symmetric 
droplet picture holds in low dimensions. Second, it is 
analytically more tractable than, say, Ising systems since 
the limit m — > oo allows for some simplifications, and 
since the replica symmetric theory is much simpler than 
the theory for broken replica symmetry. Third, it has 
been shown2*2*4 to have an unusual type of phase tran- 
sition, namely a generalized Bose-Einstein condensation 
where the spins condense in a high dimensional subspace 
in the low temperature phase. Finally, it has been stud- 
ied in a field-theoretic approach in high dimensions^Sii, 
and by computer simulations in low dimensions^. 

In this work we study the m-component spin glass on 
a Bethe lattice. The Bethe lattice allows for some exact 
analytical results using the cavity method 10 and the influ- 
ence of the finite connectivity can be analyzed explicitly. 
The cavity method has to be extended since cavity fields 
alone are not sufficient to describe the behaviour of such 
a spin system and terms of quadratic and higher order 
need to be taken into account. With this modification, 
we find a phase transition to a replica symmetric state in 
the model at a finite temperature. 

The generalized Bose-Einstein condensation is also ob- 
served on the Bethe lattice. It is characterized by an 
exponent /i which describes the scaling of the dimension- 
ality no of the ground state subspace with the number 
of spins N. For the fully connected model fj, = 2/5^; 
here we find different and connectivity-dependent values. 
In addition to /x there are a number of other zero tem- 
perature exponents: an exponent we call x describes the 
scaling of the ground state energy e(m, N) per spin and 
component with m for N = oo, the scaling of e{m,N) 
with TV for m = oo is characterized by an exponent y 
and v is the exponent of the singular part of the eigen- 
value spectrum of the inverse susceptibility matrix (for 
m = N = oo). In this paper we show that these four 
exponents are not independent but related by scaling 
laws. As noted in Ref. |j, the order of limits m — > oo 



and TV — > oo is important. Physically, taking N — > oo 
first makes most sense. Analytically, it is often more use- 
ful to take the opposite order, and this is what we will 
do in this work. Moreover, the generalized Bose-Einstein 
condensation can only be observed for the order m — ► oo 
first. The scaling relations show that even though the 
order of limits is important, it is possible to obtain infor- 
mation from one order of limits about the other. 

The paper is organized as follows. In Sec.[n]we present 
the details of our model. In Sec. lIIII we extend the cavity 
method to m-component spins, which we need in Sec. II VI 
to study the phase transition. In Sec. [V] we examine 
the generalized Bose-Einstein condensation numerically. 
The results of this are used to check the scaling relations 
derived in Sec. IVII We conclude in Sec. IVIII 



II. MODEL 

The spin glass model we analyse in this work consists 
of m-component vector spins s*j on a Bethe lattice. The 
spins have fixed length \si\ 2 = m. The Bethe lattice in 
the context of disordered systems is a random graph with 
a fixed connectivity for every spin, equal to k + 1 in this 
work. In such a graph consisting of N sites, loops are 
of order log TV. Therefore the structure of the lattice is 
locally tree like. We analyse the model in the m — > oo- 
limit for which the fully connected spin glass is replica 
symmetric^. The same is expected to hold on the Bethe 
lattice and we will restrict our work to that case. The 
generic Hamiltonian for the described model is 



H = -Y,Jir 



■ t Oj. 



(1) 



Here the exchange interactions Jij are equal to zero if 
the spins Si and Sj are not nearest neighbours on the 
Bethe lattice. If they are, Jij is drawn from a distribu- 
tion P(Jij) which is either a bimodal distribution (cor- 
responding to the ±J model) or a Gaussian distribu- 
tion. The width of these distributions J' will be cho- 
sen to scale with the number of neighbours, such that 
J' = J/yk with J fixed. This ensures easy comparison 
in the limit k — > oo with the fully-connected spin glass 
where J' = J/y/N. 
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FIG. 1: The three spins si, S2 and S3 and their back branches 
are merged onto spin s*o with coupling constants Ji , J2 and J3 
in this example with k = 3. In the next step of the iteration, 
s*o and two other branches will be merged onto the spin S4. 



case. For the time being, we write the effective Hamilto- 
nian Hi of a fc-connected spin i in general form as 

Hi = -hiSi - s^AiSi - Bi(si), (2) 

i.e. a linear term with a field hi, a quadratic term with 
a symmetric traceless matrix Ai, and the rest, contained 
in the function Bi(s) which has only terms of order s 3 
or higher. The matrix Ai may be chosen to be traceless 
without loss of generality because any nonzero trace only 
yields a constant contribution due to the constraint s$ = 
m. When k such spins are merged onto a new spin sq 
with coupling constants Ji , the partition function of this 
system is 

Z ° = J fn^5(m-sf )^ exp j-jS^-J,*,,* + 

(3) 

= f d m s 5(m- sh 2 ) f[Z l: where (4) 
^ »=i 

Zi = J cTsi 6(m - si 2 ) exp { \ n ^ +J%3o } . (5) 

Here (3 = l/ksT is the inverse temperature, as usual. 
The task is to calculate Zi and express \\ i=1 Zi in the 
form exp{— (3Hq} where Hq has the same functional form 
as in Eq. J2J. We will now proceed in two ways. The 
first is a self-consistent effective field approach where we 
use the description solely in terms of fields regardless 



III. CAVITY METHOD FOR m-COMPONENT 
SPINS 

The cavity method on a k + 1-connected Bethe lattice 
for Ising spins is explained in detail in Ref. and we 
refer the reader to this paper for a review. The main 
idea for Ising spins is that the local field at a site in 
the absence of one of its neighbours (the cavity field) is 
a quantity which can be propagated through the lattice 
iteratively. This is done by merging k branches of the 
lattice onto a new spin s , see Fig. ^ The many back 
spins in each branch may be represented by an effective 
cavity field acting on the end spins Sj. Before merging, 
the end spins Sj (i = 1, . . . , k) of these branches are each 
missing a neighbour, but after the merger they are k + 1- 
connected. The new spin sq only has k neighbours from 
the k branches that are being joined and is thus of the 
same type as the merged spins Sx,...,Sk were before the 
merger. The original spins Si, . . . , s& are then traced out, 
leaving only sq with a new effective cavity field (which 
of course depends on the original k cavity fields). This 
procedure may then be iterated, leading to a recursion re- 
lation for the cavity fields which can be used to calculate 
their distribution Q(h), see Sec. lIVI 

In the case of m-component spins this method needs to 
be modified because a simple vector field is not sufficient 
to describe the action of k old spins and their effective 
fields on a new spin. We will see below why this is the 



of the fact that it is not sufficient as mentioned above. 
The effective field will be determined by the condition 
that it has to reproduce the spin expectation value of the 
full Hamiltonian. This procedure will show that there 
is a spin glass transition at a finite temperature. The 
second way is to analyze the full problem, which is only 
possible in the vicinity of the transition. When that is 
done, it will become clear that near the transition the self- 
consistent field approach is equivalent to the first order 
approximation of the full problem. It may thus be viewed 
as the first step of a systematic approximation for the 
whole low temperature phase. 



Setting Ai = Bi(s) = 0, we have to calculate 

Zi = J d m Sl 8(m - si 2 ) exp {(3(hi + J t s )s^ . (6) 

For notational simplicity, we will abbreviate the term 
hi + JjSo by hi. We then have 

rioa+c far , ~ . 

Zi = I — d m Si exp \ j3hiS t + z(m - s 2 ) \ (7) 
= f °° +C — exp \l3 2 lj/Az + zm- — In-}. (8) 



A. Self-consistent effective field approach 
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Here we have used an integral representation of the 6- 
function and c is an arbitrary positive constant to ensure 
convergence of the Gaussian integrals over the spin com- 
ponents. The remaining integral over z may be evaluated 
in the limit m — > oo using steepest descent methods. In- 
troducing = hi/y/rn and s"q = so/^/m, the exponent 

in the integral in Eq. © is m((3 2 (fy) 2 / 4z + z - |lnz), 
apart from irrelevant constants. The saddle point equa- 
tion is therefore 

|>(/? 2 $) 2 /4* + z - \ In*) = 1 - P 2 (h 2 /^ 2 ~ ^ = 0. 
az 2 2z 

(9) 

The relevant solution of the saddle point equation is z = 



\ 1 + VI + 4/3 2 (/i<) 2 and we get 



= exp {m(/3 2 (^) 2 M ? + *° - 5 h*?) ) ■ l ' " > 



Note that according to Eq. J7J) this is nothing but the 

partition function of a single spin in a field hi . Therefore 
one can obtain the spin expectation value in a field by 

taking the derivative with respect to hi, 



1 din Zi 



Phi 
2z? 



2f3hi 



(11) 



1+ Vl + 4/3 2 ^) 2 



We will need this result below. 

We can now calculate the propagated partition func- 
tion Zq. For that purpose, the term J^so that was 

hidden in hi must be reintroduced, i.e. we have zf = 
I (l + \A + 4/? 2 (>- + J t sh') 2 ^j in the following. Ac- 
cording to Eqs. (@J and 1)10(1 . the partition function for 
spin sq is 



Zo 



I 



d m so 8(m — sq 2 ) exp 



k 

i=l 



{P 2 {h> l + J l s' a ') 2 lAz« + z°--\nz°) 



(12) 



Here it can be seen explicitly that Zq is not the partition 
function of a single spin in a field and that the cavity 
method does not close on the level of fields. The self- 
consistent approximation now consists of finding a field 
ho which generates the same spin expectation value as the 
partition function Zq. In order to calculate the latter (the 
former has already been calculated above in Eq. dJ), we 
again need to employ steepest descent methods, but this 
time not only for the auxiliary variable z from the inte- 
gral representation of the ^-function but also for all m 
spin components Sq. Here the problem arises that the 
number of integration variables, m + 1, is of the same 
order as the large parameter in the integrand, m. In 
such a situation steepest descents can not be used. The 
situation is remedied, however, by the following obser- 
vation. The spin sq only appears in the combination 
(hi + J,:So) 2 = hf + 2JihiSo + mJf, so all components of 



s*o orthogonal to the hi are projected out. These m—k or- 
thogonal components may be integrated out first as they 
are merely Gaussian integrals, leaving only k + 1 nontriv- 
ial integrals which may then be treated with the saddle 
point method. Denoting by O the orthogonal transfor- 
mation t = Osq which transforms the spin variables such 
that only the integrals over the first k components re- 
main while the rest can be carried out (the details of this 
transformation arc unimportant, as will become apparent 
below), we get (hj + J t s ) 2 = h\ + 2J l (Ohj)t + m J 2 and 

$ = \{l + V /l + 4/32($)2 + 2J i (Oh;)f' + Jf)\ with 

t ' = t/y/m. The vector t initially has m components, 
but after carrying out the integrals over the last m—k 
components only k components remain and one gets 



f dz . 
Z = — I d k t exp 

J-ioo 2tt 




{0 2 {{h' i ) 2 +2J i {Oh' i )V 



J?)/**? + *? ~ \ ^ A + *(1 - (t '?)) ~ ^ In ~ 



(13) 



The saddle point equations following from this are Since t is merely the transformed spin sq, the saddle point 

' (14) 



2(1 -(t 2 ) 

2 .r=£^. d5) 

i=l 1 
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value of t immediately gives the expectation value of the with h' — ho/ y/m. This complicated implicit equa- 



spin via (sq) = O t (padding the last to — k components 
oft with zeroes). Multiplying Eq. (fT4*|l by 2(s ) yields 

(so 



2z(s ) = 



1 - (80 



>\2 : 



(16) 



while applying O t from the left to Eq. (|15|) results in 



2z(s ) = 



/3 2 Jihi 



2z K 



(17) 



Inside z®, t may likewise be replaced by O(sq), such that 



1 



1 



4^((ft' j ) a + 2Wo , ) + J?) 
We are now in a position to present the equation which 
determines the effective field. It follows from Eq. i|ll|) 
that, if the spin sq were only subject to a field ho, 
one would have 0ha = . Combining this with 

Eqs. Jinj, JUJl and again Eq. (HJ one gets 



2/3 J* ^ 



=i 1 



l + 4/3 2 ((^) 2 +4/3J 2 



l+Vl+4^ 2 (^) 2 



(18) 



tion determines the field ho which produces the same 
spin expectation value as Zq. It corresponds to the 
Ising analogue ho — Y^%=\ u (Jii hi) from Ref . ITol where 
u(Ji 7 hi) = atanh(tanh /3Ji tanh/3/ii). In our context, 

the presence of ho on the right hand side of Eq. I|18|) 
makes the solution much more difficult. Above and close 
to the phase transition, however, an expansion of Eq. I|18|) 
in powers of hi should be valid. To first order this yields 



h 



k 



20Jihi 



(19) 



B. Analysis of the full problem 

We now turn to the full problem with the Hamiltonian 
from Eq. (j2J. As before for the self-consistent effective 

field approach, we will denote the term hi + JiSo by hi in 
the partition function Zi. The partition function Zi may 
be calculated perturbatively in the following way: 



Zi 





dz 


/ — ioo+c 


2^ 


/■200+C 


dz 


/ — -ioo+c 


2^ 


/■200+C 


dz 


' — ioo+c 


2^ 


/■200+C 


dz 


/ — ZOO + C 


2^ 



>{p(h 



d m Si exp <^ (3(hiSi + si AiSi + B 



(Si)) + z(m - sl 2 )| 



d m Si exp <j (3(hiSi - sf "(^ - AAsi) + mz ^ 



n=0 



exp < mz 



exp < mz 



J 
' z 

4"* 



.1,1 ft,, In det 



xe^{(B l ) l + -{(B 2 l ) l -(B. l ) 2 ) + 



A, 



■ hl7T 



n=0 



(20) 
(21) 

(22) 

(23) 
(24) 



The angular brackets (-)j denote the average with respect cate higher order cumulants. Now Zi from Eq. I|24(l can 



to the weight exp{[3(hiSi — Si T (^ — Ai)si)}. The expres- 
sion -| — Ai should be interpreted as a matrix expression, 
i.e. as — Ai. The dots in the last line above indi- 



again be treated by steepest descent methods in the limit 
to — ► oo , and the resulting saddle point equation is 



>-i&) T (j-->. 



-2 _ 
h' 



1 fz 

2m~p~ Tl \P ~ A 



1 d(Bj) 
to dz 



(25) 
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The quantity (Bi)i can itself be calculated by steepest de- 
scents and is given by (-Bj)j = Bi — A^j h^j + 

0(1). Note that the Hamiltonian Eq. J3J) must be ex- 
tensive in the number of degrees of freedom, which are 
the m spin components. Therefore the leading term in 
(Bi)i is of order m. The higher order cumulants can be 
calculated similarly. 

For simplicity, we will only consider the case Bi(si) = 
and hi and Ai small in the following, which is valid in the 



vicinity of the phase transition. However, an extension to 
higher orders is in principle straightforward. For small 
hi and Ai Eq. I|25|l can be solved perturbatively. For 
bookkeeping purposes a factor of e will be attached to 
hi and a factor of e 2 to A;. It will become clear below 
that this ansatz is consistent. A perturbative solution of 
Eq. I|25|l to second order in e yields 



exp J e ^ o + _^ e V/A^ - —-^1^==^^ + const. + 0(A . (26) 



Here zf° = \ (l + ^/T+IfPlf\ is the limit of zf as 

hi — » 0. From Eq. (0J Zq can be assembled and, setting 
e = 1 . the new ho and Aq can be read off. They are given 
by 

i— 1 1 

A =y -, ' A, , ' — hi ® hi ) . 

(28) 

The symbol ® denotes the tensor product. Eq. i|26|) 
shows that the new h and A are of order e resp. e 2 
such that the ansatz made above is indeed consistent. 

Comparison of Eqs. i|19f) and (|27H shows that the self- 
consistent effective field approximation and the solution 
of the full problem are identical to first order in hi. 

IV. CALCULATION OF THE PHASE 
TRANSITION 

So far we have described how to iterate on the Bethe 
lattice with a given realization of the disorder. In order 
to describe the typical behaviour, we need to average 
over the disorder. For simplicity, we will focus on the 
self-consistent effective field approach first. 

A. Phase transition at the self-consistent effective 
field level 

Instead of having a certain cavity field at each site, we 
now need to find the distribution Q(h) of cavity fields. 
Under the assumption that for m — oo the spin glass 
is replica symmetric, there is only one such distribution 



which can be found by solving the functional fixpoint 
equation 

Q(h) = Ej J (jJ^hiQih^ 8(h-h ({Ji},{hi})), 

(29) 

where h ({Ji}, {hi}) is the solution of Eq. (|T%|l . This 
equation reflects the fact that on average, all sites are 
identical such that the new distribution of fields on the 
left hand side is the same as that of the k merged spins 
under the integral on the right hand side. The symbol 
Ej stands for the average over the coupling constants 
Ji. Note that Eq. I|29|l has the trivial solution Q(h) = 
5(h). In the paramagnetic high temperature phase the 
trivial solution is expected to be stable while below the 
transition temperature it becomes unstable. This will be 
explicitly verified in the following. 

The phase transition takes place when, coming from 
high temperatures, the true local field becomes nonzero 
for the first time because it is linked to the Edwards- 
Anderson order parameter 5ea = (s) 2 via Eq. The 
true local field h t is the field which arises when joining 
k + 1 spins at a site, and its distribution is therefore 

Q t (h t ) =Ejj (nS d m h z Q(hi)) S(ti t - h ({Ji}, {hi})) 

where ho is the solution of Eq. I|18f) but with the sum 
extending to k+1 instead of k. It is thus closely related to 
the distribution Q(h) of the cavity fields, and its variance 
is zero if and only if the variance of Q(h) is zero. 

We analyze the variance of the distribution Q(h) by 
regarding Eq. (|29|1 as an iterative prescription. By mak- 
ing a Gaussian ansatz for Q(h) with width e< 1 on the 
right hand side, we calculate the width of the resulting 
distribution on the left hand side. If the new width e' 
is smaller than e, the iteration converges to 5(h), if it 
is larger, the trivial solution is unstable. Using the ap- 
proximative Eq. i|19|) . which is valid for small hi, the new 



6 



width is 

e' 2 = Af3 2 e 2 k J d.h P(Ji)^j. (30) 

Equating e and e', the critical temperature can easily be 
evaluated for the ± J- model and yields 

T c = J(l-\). (31) 

In the limit k — > oo, this agrees with the critical tem- 
perature of the fully connected model T c — Ji. For the 
Gaussian distribution of bonds, the inverse critical tem- 
perature is given by the solution of the integral equation 

1 = A(3 2 / -= 32 

J V2^J (1 + ^/l + 4(3 2 x 2 /k) 2 

which also yields T c = J in the limit k — > oo. 



B. Phase transition for the full model 

For the full problem, a similar equation as Eq. 129|) 
can be written down for the distribution function of the 
matrices A. However, while the new field ho depends only 
on the old fields hi, the matrix Aq not only depends on 
the old matrices Ai but also on the old fields, see Eq. (J2BJ. 
Therefore the question is whether the distribution of A 
becomes nontrivial at a higher temperature than Q(h). 
If it did, there would be another phase transition at that 
temperature but if it didn't, A would be slaved to the h 
and aquire a nonzero variance at the same temperature 
as Q(h). Since the calculation is essentially equal to the 
one shown in the preceding paragraph, we will merely 
quote the result that the transition temperature of A (in 
the absence of any fields) for the ± J- model would be 
T A = -^=(1 — ^Tjr) which is lower than that of the fields 
for all k, such that there is a common phase transition at 
temperature T c . This is of course physically sensible and 
agrees with the results from the fully connected model. 
The same is expected to hold for all higher orders if the 
full model is extended to include those orders. 



indistinguishable from a normal Bose-Einstein condensa- 
tion for which n = 1. We therefore proceed to analyze 
the ground state of our model numerically for finite N. 

In the ground state every spin is aligned to the direc- 
tion of its local field, 

N 

H{Si — ^ JijSj. (^^) 
j=l 

Here Hi is the modulus of the local field acting on the 
spin Si. This equation can be transformed into the eigen- 
value equation X)j=i {^Aj — Jij) s" = 0. This eigen- 
value equation for the inverse susceptibility matrix A^ = 
Hidij — Jij shows that there is at least one null eigenvalue 
of A. O n th e other hand, it was shown that A can have at 
most V2N linearly independent null eigenvectors^. Be- 
cause the dimension of the null eigenspace corresponds 
to the dimension of the space spanned by the spins 3 , 
the spins condense into a no-dimensional subspace of the 
m-dimensional space they live in with < no < y2N. 
This phenomenon was named generalized Bose-Einstein- 
condensation. The dimension no can be found by iter- 
ating Eq. i|33|) until convergence, and then diagonalizing 
the matrix A. For the fully connected model it was found 
that n ~ N» with u = 2/A 

We have repeated this computation for m-component 
spins on the Bethe lattice with k + 1 = 4, 6, 8. The accu- 
racy of reaching the groundstate was 10 -7 for the angle 
between the old and the new direction of every spin Si . A 
list of system sizes and numbers of samples used can be 
found in Table [I] An example for the eigenvalue density 
of the matrix A and the results for Ej(no) are shown in 
Figures |21 and Table [H] shows the results for the ex- 
ponent fj,. The comparison with the results for the fully 
connected models and for 2c?- and 3d-models4 puts the 
spin glass on the Bethe lattice in between those other 
models. 

TABLE I: System sizes JV and corresponding numbers of sam- 
ples D. 



N 


50 


70 


100 


150 


200 


300 


500 


1000 3000 


D 


200 


200 


150 


120 


100 


80 


50 


20 6 



V. GENERALIZED BOSE-EINSTEIN 
CONDENSATION 

So far we have analyzed the phase transition of the 
model in the thermodynamic limit N — > oo which is 
implicit in the cavity method. The generalized Bose- 
Einstein condensation, which was found for the fully 
connected models and in two- and three-dimensional 
lattices^, can however only be detected for finite N. The 
reason for this is that the dimension no of the ground 
state is of order iV M with < [i < 1, so the fraction 
no/N vanishes as N — > oo and the situation becomes 



VI. SCALING RELATIONS 

In Ref. an argument was given to explain the value 
of /i = 2/5 for no ~ N^. This value of n arises 
as the result of a competition of two energy contribu- 
tions to the ground state energy per spin and compo- 
nent e(m,N) — jMz- The first contribution to e is the 
ground state energy as a function of spin components 
in the limit N — > oo with m large but fixed. It can be 
deducecU*ii that this contribution is eo + jm~ x + 0(m~ 2 ) 
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TABLE II: Results for the exponents /i, v, x and y and the ground state energy per spin and component for different connec- 
tivities k + 1. The exponent x was calculated from fi and v using Eq. 13511 . y was calculated from Eq. 11371 . For comparison, 
the known values for the fully connected, two and three dimensional models are also shown. 



k+1 


H 


V 


X 


V 


eo 


CO 


2/5 


1/2 


1 


2/5 


-I 


8 


0.377 (± 0.003) 


0.234(±0.0003 ) 


1.34(±0.002) 


0.505 (± 0.004) 


-0.83(±0.0016) 


6 


0.362 (± 0.004) 


0.217(±0.0003 ) 


I.45(±0.002 ) 


0.525 (± 0.006) 


-0.80(±0.0015) 


4 


0.352 (± 0.005) 


0.168(±0.0013 ) 


1.57(±0.0013 ) 


0.553 (± 0.008) 


-0.73(±0.0018) 


d = 2 4 


0.33 










d = 2 4 


0.29 
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FIG. 2: Eigenvalue density for N = 300 and k + 1 = 8. The 
dashed line is proportional to (AA + A)" to illustrate that the 
initial part of p(A) scales with the exponent v. 



with eo = — 1 and x = 1. In the groundstate, no plays 
the role of an effective number of spin components as 
the spins condense into a no-dimensional subspace, thus 



this contribution is 



The sec- 



eo + |n tor m > n . 
ond contribution consists of addititonal energy costs re- 
sulting from forcing the N spins into an no-dimensional 
subspace. This contribution is important because we are 
considering the limit m —> oo before N — > oo and no is 
not fixed but proportional to N^. The energy necessary 
for this can be estimated by the amount the eigenvalue 
spectrum p(X) of A needs to be shifted in order to push no 
eigenvalues to zero, see Fig. 0J Assuming that the spec- 
trum goes as p(A) ~ X u for small A (such a behaviour was 
predicted for finite dimensional systems^), the first no 
eigenvalues occupy the space from to AA ~ , 
so the shift in energy will be of the same order. Combin- 
ing the two terms we get 



e = eo + Cl n - + c 2 (^) 1 /(^) 



(34) 



FIG. 3: Number of null eigenvalues for k + 1 = 8 as a function 
of system size N and the fit-function f(N) ~ TV 0,377 . 




P(« 



"Downward shift" 



FIG. 4: Downward shift of the semicircular eigenvalue density 
of the fully connected m-component spin glass: no of the 
smallest eigenvalues (left, shaded) become null eigenvalues 
(deltapeak, right) 



with two constants c\ and C2- Minimizing this relation 
with respect to no yields the scaling law 



1 



/' = 



x{v + \) + Y 



(35) 



Inserting the values x = 1 and v = 1/2 for the fully 
connected model, we get u — 2/5 as needed. The 
way we have presented the argument here, however, is 
more general and applies to the Bethe lattice and finite- 
dimensional lattices as well. 

For the Bethe lattice we do not know the value of x but 
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FIG. 5: The integrated eigenvalue density F(A) for N = 3000 
and k + 1 = 8. The straight line has slope v + 1 = 1.234. 



the scaling relation Eq. (|^5|l may be used to calculate it 
from /z, which can be measured as described above, and v, 
which may be obtained numerically from the eigenvalue 
spectrum p(X). For the latter, it must be kept in mind 
that we measure the spectra after the downward shift, 
i.e. the spectrum goes as p(A) ~ (A + AA)" (excluding 
the null eigenvalues), see Fig. [21 The numerical fit was 
done for N = 3000 where AA is smallest and we used 
the integrated eigenvalue density T(A) which allows for 
higher precision since there is no need for binning. An 
example of this fit is shown in Fig. [5] for k + 1 = 8. The 
results for the exponent v are shown in Table ITT1 

In addition to Eq. (|35|l we can obtain an entirely new 
scaling relation by making a scaling ansatz for the ground 
state energy per spin and component e(m, N) in the fol- 
lowing form, 



;(m, N) - e = m _a! F(mi\r'*) 



(36) 



with a scaling function F(z). When we let N — + oo while 
keeping m fixed, this yields e — eo ~ mT x as required 
by the definition of x above (provided F(0) ^ 0). On 
the other hand, when m > uq ~ 2V*, we know that 
the ground state energy becomes independent of m, so 
F(z) — const, x z x for z > no2V _/i (note that this really 
is an equality, not just an approximation). From this 
scaling ansatz it follows that e(oo, N) — eo ~ N~^ x , i.e. 
the finite size scaling of the ground state energy with the 
number of spin components m taken off to infinity before 
N scales with an exponent 



y = fxx 



1-M 
l + v' 



(37) 



The second half of this equation stems from eliminating 
x using Eq. In Fig. we show that this scaling law 
is, within numerical precision, fulfilled for k + 1 = 8. 
The same goes for the other values of k + 1 we tested 
(data not shown). Unfortunately, the data is too noisy 
to deduce the exponent y directly from it. For the fully 



FIG. 6: Ground state energy per spin and component for 
k + 1 = 8 as a function of N~ v with y calculated from the 
scaling relation Eq. I|37^ . The straight line extrapolates to the 
ground state energy eo in the limit N — > oo. 
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FIG. 7: Ground state energy per spin and component for the 
fully connected model. The straight line has slope —2/5. The 
error bars are smaller than the point size. 



connected model, Eq. l|3T|) predicts y = fi = 2/5. Fig.0 
shows that this is indeed observed. In addition to the 
exponent y, the limiting ground state energy eo can be 
obtained by extrapolating the data in Fig. [S] to N = oo. 
The results are listed in Table [HJ As k increases, eo 
converges towards —1, the value for the fully connected 
model. 



VII. CONCLUSION 

In this work we have analyzed the m-component spin 
glass on the Bethe lattice. By extending the cavity 
method, we have shown that there is a phase transi- 
tion to a replica symmetric state at a finite temperature 
for all connectivities k + 1 > 2. The phase transition 
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is a generalized Bose-Einstein condensation as for the 
fully connected modeli We have examined four differ- 
ent zero temperature exponents: fi, which describes the 
scaling of the dimension of the ground state with N (for 
?7i = oo), the exponent v of the singular part of the spec- 
trum of the inverse susceptibility matrix, x, the exponent 
which determines the scaling of the ground state energy 
per spin and component with m (for N — oo), and y, 
the scaling exponent of the same quantity with iV (for 
?7i = oo). These exponents are connected via two scaling 
laws, Eqs. (|33|) and i|37|). As we never used the partic- 
ular structure of the Bethe lattice in the derivation, we 
expect the scaling laws to hold in finite dimensions as 
well. Additionally, the exponents should be the same in 
the whole low temperature phase. The exponent u, for 
instance, should be the same at finite temperature be- 
cause it was argued in^ that the eigenvalue spectrum of 
the inverse susceptibility matrix becomes stuck as soon 
as the temperature goes below the critical temperature. 
This is also what happens in the spherical modeli. At any 
rate, the scaling laws may be used to obtain the exponent 
x, which is a quantity of interest because it extrapolates 
towards finite m for a system in the thermodynamic limit 
N = oo, from the exponents fi, y and v which are defined 
for m = oo, a case which is often computationally and 
analytically easier to handle. 



Our exponent y is analogous to the shift exponent 
1 — <d s /d as discussed in Ref. for Ising spin glasses. 
There it was found that l — Q s /d= 2/3 for all mean-field 
models considered including the Bethe lattice with con- 
nectivities 3,4,6 and 10. For the fully connected model, 
this was also confirmed by Billoire 1 ^. This is clearly not 
the case for the m = oo model, see Table ITTI For the fully 
connected model, we found y — 2/5 and increasing values 
for smaller connectivities. The origin of this difference in 
behaviour is not clear to us. 

In this work we have presented the leading contribution 
as m — ► oo. It should be pointed out, however, that our 
modification of the cavity method can easily be extended 
to obtain higher order corrections in 1/m. Work along 
these lines is in progress. 
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